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Table 3. Values of k and corresponding combinations and 
values of r 



k 

43 

67 
88 
123 
148 
173 
187 
188 
197 


Combination 
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(14-2-*) (14>3-«)L(s,x) 
(1+2-0 (H-3-0 (l+5-)£(*,x) 
(l+3-0£fox) 
(1+2-0 (l+5-0£(«,x) 
(1+3-0 (i+5-0(i+7-0£(*.x) 
(1+2-0 (1+3-0 (l+5-0£(*,x)- 
(1+2-0 (1+3-0 (l+5-0£(*,x) 
(1+3-0 i(*,x) 
(l+2-)2i(«,x) 


2 
3 
4 
2 
3 
4 
3 
3 
2 
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Forced Oscillations in Nonlinear Systems ' 

By Mary L. Cartwright 

This paper shows how the approximate form of the solutions of a certain nonlinear 
differential equation occurring in radio work may be obtained from certain general results 
and gives the proof of the general results in detail. The proof of the general statement 
depends on a type of method that can be applied with minor modifications to any equation 
of the type 

x+ k f(x) x+g(x) =kp(t), 

where p(t) has period 2tt/X, and p(t)dt is bounded for all t, f(x)>lfov \x\ >a, and 
g(x)/x> 1 for \x\ >a. 



For some years Professor J. E. Littlewood and I 
have been working on nonlinear differential equa- 
tions 2 of a type which occur in radio work and 
elsewhere. One of the most interesting of these 
equations is 



x=Jc(l — £ 2 )x+z — bk\ cos (Xt+a), 



(i) 



especially for k large and 0<6<2/3. Our attention 
was drawn to it by a remark of van der Pol, 3 which 

1 This paper contains material presented in lecture form to the staff of the 
Institute for Numerical Analysis of the National Bureau of Standards on January 
28, 1949. Miss Cartwright was a consultant at the IN A at the time this lecture 
was delivered. 

2 See M. L. Cartwright and J. E. Littlewood, J. London Math. Soc. 20, 180- 
189 (1945), and Ann. Math. 48, 472-494 (1947); also M. L. Cartwright, J. Inst. 
Elec. Eng. (Radio Section) 95 (III), 88-96 (1948, and Proc. Cambridge Phil. Soc. 
45,495 (1949). 

3 B. van der Pol, Proc. Inst. Radio Eng. 22, 1051-1086 (1934). 



suggested that it corresponded to a physical system 
investigated by him and van der Mark. 4 For 
certain values of the parameters the physical system 
had two possible stable oscillations, one of period 
Anw/\ and one of period (2n + l)27r/X. As a matter 
of fact in the case of (1), owing to the strictly 
symmetrical nonlinear function 1—x 2 , the period 
4n7r/X does not occur, but for certain values of b there 
are two stable oscillations of periods (2n ±;l)2w/\. 

It would take too long to give a complete proof 
of this statement here, but I propose to show how 
the approximate form of the solutions may be 
obtained from certain very general results, and give 
the proof of the general results in detail. The proof 



i B. van der Pol and J. van der Mark, Nature 120, 353-364 (1927) 
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of the general statement depends on the type of 
method which we use throughout. It can be 
applied with minor modifications to any equation 
of the form 



x+kf(x)x + g(x) = kp(t), 



(2) 



where p(t) has period 2?r/X, and I p(t)dt is bounded 

for all t, f(x)>l for |z|>a, and g(x)/x>l for \x\>a, 
provided of course that f, g, p satisfy the usual con- 
ditions required for the existence and uniqueness of 
solutions, and are either independent of k or satisfy 
certain simple inequalities independent of k. 

We assume throughout that k> 1, and usually 
that it is large. We first observe that the a in (1) 
is merely inserted for convenience. It enables us to 
choose the origin of t as we wish. We shall write 
x ,x for the values of x and x at / = 0. 

There are two main weapons of attack besides the 
equation (1) itself. The first is the integrated equation 



\— x +k(^— x-y+Xoj+ j vdt 
~bk [sin (\t-{-a) — sin a]; 



(3) 



the second is the energy equation which is obtained 
by multiplying (1) by 2x and integrating, 

x 2 -x 2 + 2k (' (x 2 -l)x 2 dt + x 2 -x 2 

= 2bk\ f x cos (U+a)dt (4) 

In some ways the integrated equation is more 
fundamental than (1), and that is why the X is 
inserted on the righthand side of (1). The energy 
equation may be rewritten with x=y, x =y in the 
form 

x 2 + y 2 -(xl+yl) 

= 2k\ f'&Xscos (\t + a)-(x 2 -l)x 2 \dt. (40 

The righthand side of (4') will be dominated by the 
term x 2 x 2 when x is large unless x is then very small, 
and it seems improbable that x can be small for 
most of the time that x is large. This suggests that 
x 2 -\-y 2 decreases rapidly over any arc for which x is 
large, and therefore that x is bounded for sufficiently 
larger. The general result which I propose to prove 
later is the following: 

Theorem 1: If x=x(t,x ,x ) is any solution of (1) 
for which x=x , x=x when t=0, then 



\x\<B,\x\<Bk, 



(5) 



where B is a constant independent ojk and t } for k>l, 

v y^lo \Xo) Xq) • 



Assuming for the moment the result of theorem 1 , 
we may argue as follows: eq 3 can be rearranged in 
the form 



r 3 f 1 P 

j~x=b sin (\t + a) + C-~- k \ 



xdt 



(30 



where C is a constant depending on x ,x ,b and a, 
and in virtue of (5) the last term is 0(1 /k) for large 
values of k and 0<t<27r/\. Hence there are two 
extreme possibilities: either x is comparable with k 
in magnitude, or x is given approximately by the 
equation x = X, where 

F(X)=^-X=b sin (\t + a) + C. (6) 

o 

Both these possibilities may occur for arbitrarily 
large values of t y and also of course there are transi- 
tions from one to the other. For fixed t and b and 
(7 such that ( — 2/3)<6 sin (X/ + o:) + r r <2/3, eq 6 has 
three roots, 

X 1 < -1<X 2 <1<X 3 =X&,C), 

as may be seen from figure 1. If b>2/3, 
b sin (XZ+a) + (7runs outside the interval [—2/3, 2/3] 
for some / in each period ; we suppose that 0<6<2/3 
because this is the most interesting range. It may 
be observed from figure 1 that if X x and Xz increase 
with t, X 2 decreases and vice-versa, so the solutions 
over 0<t<2ir/\ are approximately as shown iu 
figure 2, provided that \x\<B. 

If we integrate over a second period 27r/X</<47r/ 
X, we have (3') with C f in place of C where 

c '= c ~kL xdt=c -kk, ^v/+o( F ) 

8=1,2,3. 

Hence if \x\ <B, the solution is given approximately 
by X s (t,C') over the next period. The difference 
C' — C is actually of the same order as x/k, which 
occurs in the error term in (3'), but it can be shown 
that in the circumstances considered x is of the form 
x s (i,C) + 0(k~ iA ). Putting this a,ndx=X t (t,C)+0(k- 1 ) 
in (3'), and remembering that X s and x s are periodic, 
we see than the nonperiodic error is 



l£xs,odt+oQi) 




Figure 1. 
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-2 




Figure 2. 



and so the second wave will certainly be lower than 
the first in the case of X 3 by an amount of order 1/k, 
and the solution will perform descending waves as 
shown in figure 3. The second wave will be higher 
than the first in the case of X lf and higher or lower 
for X 2 , unless 



Jo 



X 2 (t 1 G)dt=0 



(so that C' = C approximately). This can only 
occur for C=0, and then we have a solution with 
period 2ttJ\ approximately. There is in fact a 
strictly periodic solution with period 2ir/\ which is, 
for reasons which we shall give presently, unstable. 

Returning to the case in which \x\ is comparable 
with k in magnitude, we suppose first that £>5fc, 
where 5>0 for a time of length D/k. Then x changes 
by at least Dd, and this is excluded by (5) if D8>2B. 
Similarly x is not less than —bk for an interval of 
time greater than DJk, for if it is x runs outside the 
strip \x\<CB. 

As regards the transitions, if |l— # 2 |>5>0 and \x\ 
is large, the second term in (1) is much more im- 
portant than x or b k\ cos(Xf+«), and so x has 
the same sign as x in |x|<^l — 8 and the opposite if 
|x|>l + 5, (provided that \x\ is large and k is large). 
Hence any small but significant deviation from X x 
or X z is at once corrected, whereas a small but signi- 
cant deviation from X 2 causes |x| to increase rapidly. 
In the latter case \x\ cannot decrease again until x ap- 
proaches either the Xi or X z curve which corresponds. 

Finally near x= ±1 the term bk\ cos (\t-\-a) dom- 
inates except when \t+ a is near nir+ 1/2 w. The dips 
of X 3 towards x=l and crests of X x toward x— — l 
do in fact occur in these intervals, and a complete 
theory of the behavior of solutions near these points 
is very complicated. However, in spite of that we 



can now form a fairly good general picture of the 
behavior of solutions. They settle above x=l into 
a long descent in waves of the form x=X z (t,C) ap- 
proximately, dipping a little lower each time until 
they reach the neighborhood of x=l. There they 
have three alternatives: (1) another wave x=X* 
just above 35=1, (2) an unstable wave x=X 2 just 
below 3?= 1, or (3) a rush down to the corresponding 
Xi which, as may be seen from figure 1, is near 
x=—2. The stable oscillations naturally follow 
the first or third alternative, and have period 
(2n-fl)27r/X or (2n—l)2ir/\ according as they have 
n+1/2 or Ti—1/2 waves above x=l. Separating 
them there are many types of unstable motion fol- 
lowing an X 2 curve which begins and ends on x = 1 for 
part or the whole of the way (or a similar X 2 curve 
after a long ascent near x= — l). From an X 2 curve 
they may pull up sharply to X x or rush down to X z 
at any stage. All types starting near x=l finally 
plunge down to the X x curve, and perform a cor- 
responding long ascent to x— — 1, the whole pheno- 
menon for x<^0 corresponding to that for x>0. All 
this can be rigorously established, and moreover the 
relative positions of solutions approximating to 
X z {t,Ci) and X z (t,C 2 ) remain the same throughout a 
long descent unless C x and C 2 differ by something 
which is extremely small for large k, such as e~ k . 

We now return to the proof of theorem 1. The 
preceding analysis depended to a large extent on the 
fact that 1— x 2 changes sign twice; in the work which 
follows the significant point is that the coefficient 
of x is positive for large x and that the function g(x) 
in (2) has the sign of x. Physicists may consider 
it intuitive that a system with a restoring force and 
positive damping for large x should have bounded 
solutions, and I hope to show that this is also intui- 
tive by mathematical commonsense. Incidentally 
the lemmas give a good deal of other information 
about the solutions. I shall as usual refer to t as 
the time, and say that t x is before or after t 2 according 
as ti<^t 2 or ti^>t 2 . The constants B are independent 
of k and t, and are not necessarily the same in each 
place unless a suffix is attached, and, as we said 
earlier fc>l. 

We first ensure that all solutions eventually come 
fairly near x=0. 

Lemma 1. A solution of (1) cannot have absolute 
value greater than 3 y * for all large t. 

.Suppose that this is not so, and that x> 3 H fort>ti. 
Integrating from t\ where x=x u x=x lj we have 

/x z x 3 \ C l 
%— Xi+kl — — x— -^+Xi J+ xdt=bk[sm(\t-{-a)—sina], 

and so 



0< f* xdt=-x + kx(l-^\+0(l) 



3 1A (t-t 



as t^ °° , where the constant implied depends on k. 
The lefthand side tends to oo ? and so x-^— o° ? but 
this implies that x-->— », which contradicts x>3 1A . 
Hence x<S V2 for arbitrarily large t } and similarly 
x> — 3 H for some largre t. 
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We next establish that as long as \x\ is not too large 
the value of \x\ cannot increase too much. 
Lemma 2. // \x\<3 H on an arc PQ, then 

\x Q \<\x P \+B 1 k ) 

and, if the arc lasts a time longer than 4/&<4,J£ Q | <2. 
More generally if \x\<b on an arc PQ, then 

\^Q\<\^p\+Bi(b )k. 

If |#|<3 H and x^>k on an arc, it can only last a 
time at most 2.3^/&<4/fr. For if it lasts longer, the 
solution travels a distance more than 2.3^ and there- 
fore cannot remain in the strip. The same is true if 
£< — k on an arc m\x\<S y \ 

Let Pi be the last point before Q at which \x\<k, 
or P itself, whichever is the latest. Suppose first 
that rr Pi >0, then ±>0 on the arc P\Q. For, if not, 
x=0 at some point of P Y Q and P x is not the last point 
at which \x\<k. By the integrated equation, since 

x Q — x Pi <B x k — xdtKBJc, 

J h 



Kffc+A. 

When x = h = 3*+hi J x=0. Integrating from the 
point Q to the point H at which x = h, we have 



Q=Xq — k H (x 2 —l)dx 



r 



xdt+bklsiniXt + a) — sin a]. 



Since x 2 -l>2 on QB, 

0<x Q -2kh l +Bk 



Suppose that x Q y>k, and that R is on x=3* so 
that x Q =x R =3 y \ On QB it is easy to see that 



-r 



x<^x Q -\-Bk — xdt 

<Bx Q -3 y H 

because x Q y>k. Also 



= x R —x Q = xdt<CBx Q t — — t 2 . 



and so t<CB d ±Q. 

The reduction in energy over any sufficiently high 
arc is established in lemma 5 which gives effect to 
the remarks we made about (4'). It is much the 
most difficult part of the proof. 

Lemma 5. If QB is an arc above x=3* beginning 
at Q on x=3 y * and ending at B on x=---X A , for a given 
B x >\ there exists a B A >B X > 1 such that if x Q >B 4 k, 



where B x >\ depends on b, which gives the result. 
Similarly if x Pi <^0, we have a^ — x Q <Bik. If the time 
is greater than 4k,\x Pl \<k, and we have the second 
form of the result. The result for \x\<b follows by L e t 
the same method. 

The next two lemmas show that the height of an 
arc QB outside the strip |a?|<3*, and the time taken Then 
to describe it are bounded by numbers depending on 
±q) in other words the velocity with which it emerges 
from the strip. 

Lemma 3. If the arc QB lies above x=3 H and 
begins at Q on x = 3 V2 , the greatest height h satisfies 



xKxI — ABMq. 
The energy eq 4 for the arc QB is 

x 2 R -x 2 Q =-2k l\x 2 -l)x 2 dt + 2bk\ f x sin (\t + a)dt. 

Jo Jo 



J=k I ' (x 2 -l)x 2 dt>k Cx 2 dt. 

Jo Jo 



(8) 



4-4<-2J-f 2bk\ f \x\dt 
<-2J+BM H (l t \i\ 2 dty 
by Cauchy's inequality. So by lemma 4 and (8) 

4-4<-2J+.B 5 fc*s$<7* 

<-J<-BlH Q , 

provided that J>B? 5 kx Q , where B 5 depends on B 3 and 
x Q >k. We may obviously choose 5?>4J?i, and we 
have the result except for the case in which 
J<Bikx Q . 

Suppose now that J<^B 2 5 kx Q , and integrate from 

Q until x=- x Q or x = 2B\-\-?> v % whichever comes 



and the result follows for hi and so for h. 

Lemma 4. If QB is an arc above x=8 H beginning this arc and 
at Q on x—3 H j then the time t taken to describe QB is 
less than B z x Q , provided that x Q y>k. 



first. Then x> — ± o >0. 



Hence x is increasing on 



f sin Qit + a)xdt= \ sin (\t + a)dx< | dx<2B\. 
Jo Jz y * Jz y * 
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Hence 



x 2 -x 2 Q =-2k f\x 2 -l)x 2 dt-kb\ f\m(\t+a)dx-B 
Jo Jo 



>-2J-kBB\-B 
>-2kB 2 5 x Q -kB>-^x 2 Q 



provided that i >&Z? 4 . It follows that x 2 > -x% and 

therefore £>-£ Q , so that x reaches 2B 2 5 -\-3 H first. 
Now since x 2 — 1>2 on the arc, 

r 255+3 H 
J>2k \x\dx>kx Q '2Bi, 

Js y * 

and we have a contradiction. Hence the result of 
the lemma is true. 

It remains to combine the results of lemmas 2 
and 5. 

Lemma 6. If Q is a point on x=3 1A such that 
x Q ^>B^k^>0 } then the solution returns to the strip 
\x\<3 y * at R and emerges again at S with \x s \<C\x Q \ 
— Bik. 

By lemma 2 with R and S in place of P and Q 

\x s \ 2 <\x R \ 2 +2B l k\x R \+B 2 k 2 . 




Using lemma 5 in the form stated and also in the 
weaker form |#22|<|xq|, we have 

l^| 2 <l^| 2 ~45iA:!^| + 25 1 /:fe|+jB^ 

= {\x Q \-B x k)\ 

From which the result follows. For x Q ^>BJc> B x k. 
We now have the result stated in theorem 1 . For 
by lemma 1 the solution must enter the strip 
|#|<3 H , and by the second part of lemma 2 \x\<B x k y 
if it stays there. If on the other hand it emerges 
at Q with \x Q \<BJc, the height (or depth) of the 
subsequent arc outside the strip is less than B A -\-B 2 . 
For by symmetry all the lemmas for x— >3 1/a have 
strictly corresponding forms for x< — 3 H . But by 
lemma 6 if \x Q \y>B 4 k, the solution emerges the next 
time at S with |x s |<ji Q j— B x k, so that it must either 
stay in the strip or emerge eventually with \x Q \< BJc y 
and then by the second part of lemma 2, \x\ < B±-\-B 2 
for all subsequent t. 
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Figure 5. 
Los Angeles, December 16, 1949. 
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